Unbiased in silico design of pH-sensitive tetrapeptides.

We used coarse-grain molecular dynamics simulations to screen all possible histidine-bearing tetrapeptide sequences, finding novel peptide sequences with pH-tunable assembly properties. These tetrapeptides could be used for various biological applications, such as triggered delivery of bioactive molecules.

Molecular determinants of TRAF6 binding specificity suggest that native interaction partners are not optimized for affinity.

TRAF6 is an adaptor protein involved in signaling pathways that are essential for development and the immune system. It participates in many protein-protein interactions, some of which are mediated by the C-terminal MATH domain, which binds to short peptide segments containing the motif PxExx[FYWHDE], where x is any amino acid. Blocking MATH domain interactions is associated with favorable effects in various disease models. To better define TRAF6 MATH domain binding preferences, we screened a combinatorial library using bacterial cell-surface peptide display. We identified 236 of the best TRAF6-interacting peptides and a set of 1,200 peptides that match the sequence PxE but do not bind TRAF6 MATH. The peptides that were most enriched in the screen bound TRAF6 tighter than previously measured native peptides. To better understand the structural basis for TRAF6 interaction preferences, we built all-atom structural models of the MATH domain in complex with high-affinity binders and nonbinders identified in the screen. We identified favorable interactions for motif features in binders as well as negative design elements distributed across the motif that can disfavor or preclude binding. Searching the human proteome revealed that the most biologically relevant TRAF6 motif matches occupy a different sequence space from the best hits discovered in combinatorial library screening, suggesting that native interactions are not optimized for affinity. Our experimentally determined binding preferences and structural models support the design of peptide-based interaction inhibitors with higher affinities than endogenous TRAF6 ligands.

In Silico Engineering of Enzyme Access Tunnels.

Enzyme engineering is a tailoring process that allows the modification of naturally occurring enzymes to provide them with improved catalytic efficiency, stability, or specificity. By introducing partial modifications to their sequence and to their structural features, enzyme engineering can transform natural enzymes into more efficient, specific and resistant biocatalysts and render them suitable for virtually countless industrial processes. Current enzyme engineering methods mostly target the active site of the enzyme, where the catalytic reaction takes place. Nonetheless, the tunnel that often connects the surface of an enzyme with its buried active site plays a key role in the activity of the enzyme as it acts as a gatekeeper and regulates the access of the substrate to the catalytic pocket. Hence, there is an increasing interest in targeting the sequence and the structure of substrate entrance tunnels in order to fine-tune enzymatic activity, regulate substrate specificity, or control reaction promiscuity.In this chapter, we describe the use of a rational in silico design and screening method to engineer the access tunnel of a fructosyl peptide oxidase with the aim to facilitate access to its catalytic site and to expand its substrate range. Our goal is to engineer this class of enzymes in order to utilize them for the direct detection of glycated proteins in diabetes monitoring devices. The design strategy involves remodeling of the backbone structure of the enzyme , a feature that is not possible with conventional enzyme engineering techniques such as single-point mutagenesis and that is highly unlikely to occur using a directed evolution approach.The proposed strategy, which results in a significant reduction in cost and time for the experimental production and characterization of candidate enzyme variants, represents a promising approach to the expedited identification of novel and improved enzymes. Rational enzyme design aims to provide in silico strategies for the fast, accurate, and inexpensive development of biocatalysts that can meet the needs of multiple industrial sectors, thus ultimately promoting the use of green chemistry and improving the efficiency of chemical processes.

Rational backbone redesign of a fructosyl peptide oxidase to widen its active site access tunnel.

Fructosyl peptide oxidases (FPOXs) are enzymes currently used in enzymatic assays to measure the concentration of glycated hemoglobin and albumin in blood samples, which serve as biomarkers of diabetes. However, since FPOX are unable to work directly on glycated proteins, current enzymatic assays are based on a preliminary proteolytic digestion of the target proteins. Herein, to improve the speed and costs of the enzymatic assays for diabetes testing, we applied a rational design approach to engineer a novel enzyme with a wider access tunnel to the catalytic site, using a combination of Rosetta design and molecular dynamics simulations. Our final design, L3_35A, shows a significantly wider and shorter access tunnel, resulting from the deletion of five-amino acids lining the gate structures and from a total of 35 point mutations relative to the wild-type (WT) enzyme. Indeed, upon experimental testing, our engineered enzyme shows good structural stability and maintains significant activity relative to the WT.

The Anti-Amyloidogenic Action of Doxycycline: A Molecular Dynamics Study on the Interaction with Aß42.

The pathological aggregation of amyloidogenic proteins is a hallmark of many neurological diseases, including Alzheimer's disease and prion diseases. We have shown both in vitro and in vivo that doxycycline can inhibit the aggregation of Aß42 amyloid fibrils and disassemble mature amyloid fibrils. However, the molecular mechanisms of the drug's anti-amyloidogenic property are not understood. In this study, a series of molecular dynamics simulations were performed to explain the molecular mechanism of the destabilization of Aß42 fibrils by doxycycline and to compare the action of doxycycline with those of iododoxorubicin (a toxic structural homolog of tetracyclines), curcumin (known to have anti-amyloidogenic activity) and gentamicin (an antibiotic with no experimental evidence of anti-amyloidogenic properties). We found that doxycycline tightly binds the exposed hydrophobic amino acids of the Aß42 amyloid fibrils, partly leading to destabilization of the fibrillar structure. Clarifying the molecular determinants of doxycycline binding to Aß42 may help devise further strategies for structure-based drug design for Alzheimer's disease.

Molecular dynamics investigation of halogenated amyloidogenic peptides.

Besides their biomolecular relevance, amyloids, generated by the self-assembly of peptides and proteins, are highly organized structures useful for nanotechnology applications. The introduction of halogen atoms in these peptides, and thus the possible formation of halogen bonds, allows further possibilities to finely tune the amyloid nanostructure. In this work, we performed molecular dynamics simulations on different halogenated derivatives of the ß-amyloid peptide core-sequence KLVFF, by using a modified AMBER force field in which the σ-hole located on the halogen atom is modeled with a positively charged extra particle. The analysis of equilibrated structures shows good agreement with crystallographic data and experimental results, in particular concerning the formation of halogen bonds and the stability of the supramolecular structures. The modified force field described here allows describing the atomistic details contributing to peptides aggregation, with particular focus on the role of halogen bonds. This framework can potentially help the design of novel halogenated peptides with desired aggregation propensity. Graphical abstract Molecular dynamics investigation of halogenated amyloidogenic peptides.

Thermal stabilization of the deglycating enzyme Amadoriase I by rational design.

Amadoriases are a class of FAD-dependent enzymes that are found in fungi, yeast and bacteria and that are able to hydrolyze glycated amino acids, cleaving the sugar moiety from the amino acidic portion. So far, engineered Amadoriases have mostly found practical application in the measurement of the concentration of glycated albumin in blood samples. However, these engineered forms of Amadoriases show relatively low absolute activity and stability levels, which affect their conditions of use. Therefore, enzyme stabilization is desirable prior to function-altering molecular engineering. In this work, we describe a rational design strategy based on a computational screening method to evaluate a library of potentially stabilizing disulfide bonds. Our approach allowed the identification of two thermostable Amadoriase I mutants (SS03 and SS17) featuring a significantly higher T50 (55.3 °C and 60.6 °C, respectively) compared to the wild-type enzyme (52.4 °C). Moreover, SS17 shows clear hyperstabilization, with residual activity up to 95 °C, whereas the wild-type enzyme is fully inactive at 55 °C. Our computational screening method can therefore be considered as a promising approach to expedite the design of thermostable enzymes.

Review: Engineering of thermostable enzymes for industrial applications.

The catalytic properties of some selected enzymes have long been exploited to carry out efficient and cost-effective bioconversions in a multitude of research and industrial sectors, such as food, health, cosmetics, agriculture, chemistry, energy, and others. Nonetheless, for several applications, naturally occurring enzymes are not considered to be viable options owing to their limited stability in the required working conditions. Over the years, the quest for novel enzymes with actual potential for biotechnological applications has involved various complementary approaches such as mining enzyme variants from organisms living in extreme conditions (extremophiles), mimicking evolution in the laboratory to develop more stable enzyme variants, and more recently, using rational, computer-assisted enzyme engineering strategies. In this review, we provide an overview of the most relevant enzymes that are used for industrial applications and we discuss the strategies that are adopted to enhance enzyme stability and/or activity, along with some of the most relevant achievements. In all living species, many different enzymes catalyze fundamental chemical reactions with high substrate specificity and rate enhancements. Besides specificity, enzymes also possess many other favorable properties, such as, for instance, cost-effectiveness, good stability under mild pH and temperature conditions, generally low toxicity levels, and ease of termination of activity. As efficient natural biocatalysts, enzymes provide great opportunities to carry out important chemical reactions in several research and industrial settings, ranging from food to pharmaceutical, cosmetic, agricultural, and other crucial economic sectors.

Nanostructure and stability of calcitonin amyloids.

Calcitonin is a 32-amino acid thyroid hormone that can form amyloid fibrils. The structural basis of the fibril formation and stabilization is still debated and poorly understood. The reason is that NMR data strongly suggest antiparallel ß-sheet calcitonin assembly, whereas modeling studies on the short DFNKF peptide (corresponding to the sequence from Asp15 to Phe19 of human calcitonin and reported as the minimal amyloidogenic module) show that it assembles with parallel ß-sheets. In this work, we first predict the structure of human calcitonin through two complementary molecular dynamics (MD) methods, finding that human calcitonin forms an α-helix. We use extensive MD simulations to compare previously proposed calcitonin fibril structures. We find that two conformations, the parallel arrangement and one of the possible antiparallel structures (with Asp15 and Phe19 aligned), are highly stable and ordered. Nonetheless, fibrils with parallel molecules show bulky loops formed by residues 1 to 7 located on the same side, which could limit or prevent the formation of larger amyloids. We investigate fibrils formed by the DFNKF peptide by simulating different arrangements of this amyloidogenic core sequence. We show that DFNKF fibrils are highly stable when assembled in parallel ß-sheets, whereas they quickly unfold in antiparallel conformation. Our results indicate that the DFNKF peptide represents only partially the full-length calcitonin behavior. Contrary to the full-length polypeptide, in fact, the DFNKF sequence is not stable in antiparallel conformation, suggesting that the residue flanking the amyloidogenic peptide contributes to the stabilization of the experimentally observed antiparallel ß-sheet packing.

Molecular dynamics simulations provide insights into the substrate specificity of FAOX family members.

Enzymatic assays based on Fructosyl Amino Acid Oxidases (FAOX) represent a potential, rapid and economical strategy to measure glycated hemoglobin (HbA1c), which is in turn a reliable method to monitor the insurgence and the development of diabetes mellitus. However, the engineering of naturally occurring FAOX to specifically recognize fructosyl-valine (the glycated N-terminal residue of HbA1c) has been hindered by the paucity of information on the tridimensional structures and catalytic residues of the different FAOX that exist in nature, and in general on the molecular mechanisms that regulate specificity in this class of enzymes. In this study, we use molecular dynamics simulations and advanced modeling techniques to investigate five different relevant wild-type FAOX (Amadoriase I, Amadoriase II, PnFPOX, FPOX-E and N1-1-FAOD) in order to elucidate the molecular mechanisms that drive their specificity towards polar and nonpolar substrates. Specifically, we compare these five different FAOX in terms of overall folding, ligand entry tunnels, ligand binding residues and ligand binding energies. Our work will contribute to future enzyme structure modifications aimed at the rational design of novel biosensors for the monitoring of blood glucose levels.

Crystal structure of the deglycating enzyme Amadoriase I in its free form and substrate-bound complex.

Amadoriases, also known as fructosyl amine oxidases (FAOX), are enzymes that catalyze the de-glycosylation of fructosyl amino acids. As such, they are excellent candidates for the development of enzyme-based diagnostic and therapeutic tools against age- and diabetes-induced protein glycation. However, mostly because of the lack of a complete structural characterization of the different members of the family, the molecular bases of their substrate specificity have yet to be fully understood. The high resolution crystal structures of the free and the substrate-bound form of Amadoriase I shown herein allow for the identification of key structural features that account for the diverse substrate specificity shown by this class of enzymes. This is of particular importance in the context of the rather limited and partially incomplete structural information that has so far been available in the literature on the members of the FAOX family. Moreover, using molecular dynamics simulations, we describe the tunnel conformation and the free energy profile experienced by the ligand in going from bulk water to the catalytic cavity, showing the presence of four gating helices/loops, followed by an "L-shaped" narrow cavity. In summary, the tridimensional architecture of Amadoriase I presented herein provides a reference structural framework for the design of novel enzymes for diabetes monitoring and protein deglycation. Proteins 2016; 84:744-758. © 2016 Wiley Periodicals, Inc.